Rockfish Historical Harvest Estimates

This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.

Issue	Howard	Bayes
Time series	1999 - present	1977 - present
Bias in SWHS	Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests.	Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases	Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data.	Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations	Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas.	Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation	Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with.	By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Results

• Model: LB_fit_3pH

• Model Run: Gen4int_indcomp_swhsR_FULL_pHB4pars_re0d_thru2024_4e+05_2025-10-29 (Naming convention: modelname_thruYEAR_iters_rundate)

• Proportion of parmaters converged (Rhat <= 1.1): 0.9959

• See end of document for summary of unconverged parameters.

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

And black rockfish releases…

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Total Biomass Removal Estimates

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Model fit

Logbook residuals

Figure 20.- Residuals from logbook harvests.

SWHS residuals

Figure 21.- Residuals from SWHS harvests.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

Figure 23.- Mean percent of harvest by charter anglers.

When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(slope|non-pelagic & non-yellowye) For release estimates

## NULL

Weight Fits

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.3821138

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.11 while 0.89 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable	area_n	year	user	mean	median	lower	upper	Rhat	sd	skew	area
p_pelagic	11	1	1	0.458	0.458	0.186	0.723	1.102	0.140	0.000	CSEO
p_pelagic	11	1	2	0.391	0.386	0.144	0.655	1.094	0.132	0.038	CSEO
p_pelagic	11	10	1	0.455	0.454	0.196	0.716	1.099	0.133	0.008	CSEO
p_pelagic	11	10	2	0.391	0.388	0.165	0.652	1.100	0.128	0.023	CSEO
p_pelagic	11	11	1	0.456	0.455	0.200	0.709	1.101	0.133	0.008	CSEO
p_pelagic	11	11	2	0.389	0.383	0.154	0.645	1.093	0.128	0.047	CSEO
p_pelagic	11	12	1	0.459	0.459	0.193	0.722	1.097	0.134	0.000	CSEO
p_pelagic	11	12	2	0.394	0.390	0.160	0.649	1.086	0.129	0.031	CSEO
p_pelagic	11	13	1	0.459	0.458	0.207	0.715	1.112	0.132	0.008	CSEO
p_pelagic	11	13	2	0.389	0.385	0.158	0.646	1.091	0.128	0.031	CSEO
p_pelagic	11	14	1	0.472	0.476	0.225	0.724	1.109	0.131	-0.031	CSEO
p_pelagic	11	14	2	0.397	0.394	0.161	0.658	1.082	0.131	0.023	CSEO
p_pelagic	11	15	1	0.429	0.425	0.184	0.690	1.098	0.132	0.030	CSEO
p_pelagic	11	15	2	0.387	0.382	0.161	0.644	1.088	0.125	0.040	CSEO
p_pelagic	11	16	1	0.476	0.476	0.226	0.724	1.078	0.126	0.000	CSEO
p_pelagic	11	16	2	0.396	0.391	0.162	0.653	1.091	0.127	0.039	CSEO
p_pelagic	11	17	1	0.467	0.466	0.229	0.716	1.078	0.127	0.008	CSEO
p_pelagic	11	17	2	0.394	0.388	0.172	0.651	1.080	0.125	0.048	CSEO
p_pelagic	11	18	1	0.472	0.470	0.234	0.723	1.083	0.127	0.016	CSEO
p_pelagic	11	18	2	0.399	0.396	0.172	0.669	1.072	0.127	0.024	CSEO
p_pelagic	11	19	1	0.475	0.474	0.243	0.721	1.066	0.122	0.008	CSEO
p_pelagic	11	19	2	0.400	0.391	0.178	0.654	1.057	0.123	0.073	CSEO
p_pelagic	11	2	1	0.457	0.456	0.184	0.715	1.094	0.137	0.007	CSEO
p_pelagic	11	2	2	0.388	0.385	0.147	0.656	1.123	0.132	0.023	CSEO
p_pelagic	11	20	1	0.477	0.474	0.238	0.712	1.048	0.126	0.024	CSEO
p_pelagic	11	20	2	0.403	0.397	0.170	0.669	1.050	0.126	0.047	CSEO
p_pelagic	11	21	1	0.485	0.483	0.251	0.727	1.042	0.124	0.016	CSEO
p_pelagic	11	21	2	0.403	0.398	0.173	0.655	1.028	0.124	0.040	CSEO
p_pelagic	11	22	2	0.411	0.408	0.186	0.663	1.036	0.125	0.024	CSEO
p_pelagic	11	23	2	0.448	0.443	0.220	0.693	1.027	0.123	0.041	CSEO
p_pelagic	11	3	1	0.458	0.463	0.193	0.707	1.108	0.135	-0.037	CSEO
p_pelagic	11	3	2	0.391	0.382	0.151	0.665	1.092	0.132	0.068	CSEO
p_pelagic	11	4	1	0.457	0.456	0.185	0.732	1.092	0.139	0.007	CSEO
p_pelagic	11	4	2	0.389	0.384	0.156	0.656	1.107	0.130	0.039	CSEO
p_pelagic	11	5	1	0.453	0.453	0.190	0.710	1.084	0.136	0.000	CSEO
p_pelagic	11	5	2	0.388	0.384	0.147	0.653	1.085	0.132	0.030	CSEO
p_pelagic	11	6	1	0.455	0.453	0.201	0.708	1.087	0.131	0.015	CSEO
p_pelagic	11	6	2	0.392	0.386	0.145	0.652	1.092	0.130	0.046	CSEO
p_pelagic	11	7	1	0.451	0.450	0.188	0.710	1.091	0.137	0.007	CSEO
p_pelagic	11	7	2	0.393	0.389	0.152	0.652	1.119	0.132	0.030	CSEO
p_pelagic	11	8	1	0.456	0.458	0.193	0.720	1.094	0.137	-0.015	CSEO
p_pelagic	11	8	2	0.392	0.386	0.154	0.646	1.068	0.130	0.046	CSEO
p_pelagic	11	9	1	0.455	0.453	0.197	0.717	1.105	0.133	0.015	CSEO
p_pelagic	11	9	2	0.386	0.382	0.151	0.655	1.101	0.130	0.031	CSEO
p_yellow	13	15	1	0.420	0.367	0.125	0.987	1.010	0.222	0.239	NSEI
p_yellow	13	19	2	0.423	0.374	0.129	0.986	1.010	0.219	0.224	NSEI
p_yellow	13	2	1	0.423	0.372	0.131	0.986	1.011	0.219	0.232	NSEI
p_yellow	13	20	1	0.420	0.368	0.130	0.988	1.010	0.217	0.240	NSEI
p_yellow	13	22	1	0.421	0.367	0.125	0.987	1.012	0.218	0.248	NSEI
p_yellow	13	22	2	0.419	0.368	0.129	0.986	1.013	0.218	0.234	NSEI
p_yellow	13	23	1	0.387	0.333	0.106	0.984	1.011	0.223	0.242	NSEI
p_yellow	13	24	1	0.388	0.334	0.102	0.986	1.012	0.222	0.243	NSEI
p_yellow	13	25	1	0.390	0.338	0.107	0.982	1.011	0.220	0.236	NSEI
p_yellow	13	26	1	0.382	0.324	0.108	0.984	1.013	0.219	0.264	NSEI
p_yellow	13	27	1	0.373	0.317	0.105	0.984	1.015	0.216	0.259	NSEI
p_yellow	13	28	1	0.380	0.332	0.107	0.980	1.013	0.212	0.226	NSEI
p_yellow	13	28	2	0.466	0.437	0.152	0.988	1.012	0.209	0.139	NSEI
p_yellow	13	29	1	0.359	0.315	0.097	0.961	1.017	0.204	0.215	NSEI
p_yellow	13	29	2	0.457	0.426	0.154	0.974	1.011	0.201	0.154	NSEI
p_slope	14	1	2	0.071	0.053	0.010	0.238	1.011	0.063	0.287	NSEO
p_dsr	14	45	1	0.078	0.058	0.013	0.269	1.015	0.068	0.295	NSEO
p_slope	14	6	2	0.074	0.055	0.010	0.250	1.011	0.065	0.292	NSEO
p_pelagic	16	1	1	0.388	0.383	0.136	0.676	1.142	0.144	0.035	SSEO
p_pelagic	16	1	2	0.365	0.358	0.127	0.648	1.157	0.136	0.051	SSEO
p_pelagic	16	10	1	0.398	0.391	0.152	0.679	1.124	0.139	0.050	SSEO
p_pelagic	16	10	2	0.376	0.362	0.139	0.664	1.113	0.140	0.100	SSEO
p_pelagic	16	11	1	0.398	0.394	0.153	0.682	1.110	0.141	0.028	SSEO
p_pelagic	16	11	2	0.375	0.369	0.135	0.662	1.108	0.136	0.044	SSEO
p_pelagic	16	12	1	0.401	0.395	0.157	0.693	1.104	0.139	0.043	SSEO
p_pelagic	16	12	2	0.375	0.364	0.131	0.666	1.110	0.140	0.079	SSEO
p_pelagic	16	13	1	0.401	0.401	0.153	0.674	1.118	0.136	0.000	SSEO
p_pelagic	16	13	2	0.373	0.368	0.136	0.644	1.124	0.133	0.038	SSEO
p_pelagic	16	14	1	0.384	0.378	0.148	0.658	1.096	0.133	0.045	SSEO
p_pelagic	16	14	2	0.376	0.366	0.151	0.658	1.114	0.133	0.075	SSEO
p_pelagic	16	15	1	0.394	0.388	0.160	0.666	1.108	0.135	0.045	SSEO
p_pelagic	16	15	2	0.381	0.372	0.145	0.664	1.103	0.136	0.066	SSEO
p_pelagic	16	16	1	0.383	0.376	0.162	0.651	1.073	0.129	0.054	SSEO
p_pelagic	16	16	2	0.378	0.373	0.147	0.650	1.083	0.131	0.038	SSEO
p_pelagic	16	17	1	0.443	0.440	0.203	0.701	1.074	0.131	0.023	SSEO
p_pelagic	16	17	2	0.397	0.392	0.160	0.661	1.080	0.132	0.038	SSEO
p_pelagic	16	18	1	0.421	0.418	0.195	0.670	1.052	0.126	0.024	SSEO
p_pelagic	16	18	2	0.400	0.394	0.169	0.664	1.049	0.130	0.046	SSEO
p_pelagic	16	19	1	0.410	0.405	0.179	0.672	1.037	0.128	0.039	SSEO
p_pelagic	16	19	2	0.394	0.384	0.161	0.659	1.046	0.129	0.078	SSEO
p_pelagic	16	2	1	0.386	0.382	0.139	0.669	1.127	0.141	0.028	SSEO
p_pelagic	16	2	2	0.367	0.359	0.127	0.660	1.134	0.140	0.057	SSEO
p_pelagic	16	20	1	0.416	0.409	0.190	0.683	1.037	0.127	0.055	SSEO
p_pelagic	16	20	2	0.399	0.392	0.173	0.656	1.029	0.126	0.056	SSEO
p_pelagic	16	21	1	0.400	0.392	0.186	0.656	1.030	0.122	0.066	SSEO
p_pelagic	16	21	2	0.394	0.388	0.173	0.645	1.033	0.126	0.048	SSEO
p_pelagic	16	22	2	0.402	0.396	0.181	0.661	1.018	0.125	0.048	SSEO
p_pelagic	16	23	2	0.317	0.306	0.136	0.567	1.014	0.111	0.099	SSEO
p_pelagic	16	3	1	0.386	0.382	0.137	0.660	1.166	0.140	0.029	SSEO
p_pelagic	16	3	2	0.365	0.356	0.127	0.649	1.128	0.139	0.065	SSEO
p_dsr	16	31	1	0.406	0.406	0.376	0.437	1.336	0.016	0.000	SSEO
p_pelagic	16	4	1	0.389	0.382	0.138	0.664	1.159	0.140	0.050	SSEO
p_pelagic	16	4	2	0.374	0.366	0.131	0.649	1.128	0.140	0.057	SSEO
p_pelagic	16	41	1	0.794	0.794	0.786	0.801	1.110	0.004	0.000	SSEO
p_yellow	16	41	1	0.430	0.430	0.409	0.451	1.046	0.011	0.000	SSEO
p_dsr	16	41	1	0.488	0.481	0.448	0.538	3.077	0.026	0.273	SSEO
p_slope	16	41	1	0.115	0.125	0.071	0.148	4.755	0.025	-0.406	SSEO
p_pelagic	16	5	1	0.392	0.386	0.141	0.679	1.135	0.142	0.042	SSEO
p_pelagic	16	5	2	0.371	0.364	0.131	0.658	1.150	0.138	0.051	SSEO
p_pelagic	16	6	1	0.390	0.384	0.142	0.671	1.129	0.140	0.043	SSEO
p_pelagic	16	6	2	0.367	0.362	0.125	0.645	1.132	0.138	0.036	SSEO
p_pelagic	16	7	1	0.389	0.382	0.143	0.668	1.147	0.140	0.050	SSEO
p_pelagic	16	7	2	0.371	0.366	0.131	0.653	1.114	0.137	0.037	SSEO
p_pelagic	16	8	1	0.394	0.391	0.152	0.668	1.123	0.135	0.022	SSEO
p_pelagic	16	8	2	0.373	0.369	0.130	0.659	1.141	0.138	0.029	SSEO
p_pelagic	16	9	1	0.396	0.392	0.149	0.672	1.117	0.138	0.029	SSEO
p_pelagic	16	9	2	0.368	0.363	0.127	0.645	1.125	0.137	0.037	SSEO
p_pelagic	3	1	1	0.552	0.556	0.285	0.797	1.020	0.130	-0.031	PWSI
p_pelagic	3	1	2	0.188	0.175	0.064	0.386	1.014	0.084	0.155	PWSI
p_pelagic	3	10	1	0.554	0.560	0.288	0.789	1.015	0.132	-0.046	PWSI
p_pelagic	3	10	2	0.189	0.176	0.067	0.382	1.026	0.082	0.158	PWSI
p_pelagic	3	11	1	0.549	0.551	0.286	0.790	1.019	0.132	-0.015	PWSI
p_pelagic	3	11	2	0.190	0.177	0.065	0.391	1.016	0.084	0.156	PWSI
p_pelagic	3	12	1	0.556	0.558	0.302	0.788	1.026	0.128	-0.016	PWSI
p_pelagic	3	12	2	0.190	0.174	0.067	0.396	1.012	0.086	0.187	PWSI
p_pelagic	3	13	1	0.553	0.556	0.297	0.784	1.016	0.128	-0.023	PWSI
p_pelagic	3	13	2	0.190	0.177	0.066	0.392	1.012	0.084	0.156	PWSI
p_pelagic	3	14	1	0.550	0.555	0.290	0.789	1.011	0.130	-0.038	PWSI
p_pelagic	3	14	2	0.195	0.182	0.067	0.388	1.016	0.086	0.152	PWSI
p_pelagic	3	15	1	0.546	0.550	0.283	0.794	1.017	0.132	-0.030	PWSI
p_pelagic	3	15	2	0.187	0.174	0.066	0.382	1.013	0.083	0.157	PWSI
p_pelagic	3	16	1	0.547	0.551	0.289	0.777	1.010	0.129	-0.031	PWSI
p_pelagic	3	16	2	0.189	0.176	0.071	0.393	1.012	0.083	0.157	PWSI
p_pelagic	3	17	1	0.555	0.562	0.293	0.792	1.022	0.131	-0.054	PWSI
p_pelagic	3	17	2	0.188	0.175	0.065	0.375	1.019	0.082	0.159	PWSI
p_pelagic	3	18	1	0.553	0.560	0.288	0.793	1.013	0.131	-0.053	PWSI
p_pelagic	3	18	2	0.190	0.176	0.067	0.388	1.013	0.083	0.169	PWSI
p_pelagic	3	19	1	0.553	0.555	0.308	0.791	1.018	0.127	-0.016	PWSI
p_pelagic	3	2	1	0.551	0.559	0.287	0.789	1.016	0.133	-0.060	PWSI
p_pelagic	3	2	2	0.192	0.181	0.067	0.384	1.016	0.084	0.132	PWSI
p_pelagic	3	3	1	0.549	0.555	0.292	0.785	1.013	0.128	-0.047	PWSI
p_pelagic	3	3	2	0.189	0.177	0.063	0.391	1.016	0.084	0.142	PWSI
p_pelagic	3	4	1	0.548	0.549	0.291	0.790	1.026	0.132	-0.008	PWSI
p_pelagic	3	4	2	0.190	0.177	0.067	0.389	1.013	0.086	0.152	PWSI
p_pelagic	3	5	1	0.551	0.555	0.283	0.787	1.013	0.131	-0.031	PWSI
p_pelagic	3	5	2	0.189	0.176	0.067	0.387	1.015	0.083	0.157	PWSI
p_pelagic	3	6	1	0.553	0.559	0.295	0.778	1.021	0.125	-0.048	PWSI
p_pelagic	3	6	2	0.190	0.178	0.066	0.383	1.012	0.083	0.145	PWSI
p_pelagic	3	7	1	0.552	0.557	0.290	0.792	1.022	0.130	-0.039	PWSI
p_pelagic	3	7	2	0.192	0.178	0.066	0.400	1.011	0.085	0.165	PWSI
p_pelagic	3	8	1	0.548	0.551	0.286	0.790	1.020	0.131	-0.023	PWSI
p_pelagic	3	8	2	0.190	0.177	0.070	0.393	1.015	0.085	0.153	PWSI
p_pelagic	3	9	1	0.550	0.557	0.290	0.779	1.022	0.129	-0.054	PWSI
p_pelagic	3	9	2	0.190	0.177	0.066	0.393	1.016	0.085	0.153	PWSI
p_pelagic	4	14	1	0.539	0.545	0.276	0.774	1.010	0.131	-0.046	PWSO
p_pelagic	4	2	2	0.540	0.545	0.284	0.779	1.010	0.129	-0.039	PWSO
p_black	8	40	1	0.458	0.414	0.342	0.956	1.090	0.155	0.284	afognak
p_black	8	48	1	0.540	0.492	0.412	0.971	1.014	0.151	0.318	afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable	area_n	year/species	mean	median	lower	upper	Rhat	sd	skew	year	area
beta0_pH	11	1	0.118	0.144	-0.610	0.696	1.315	0.356	-0.073		CSEO
beta1_pH	11	1	3.273	3.212	2.414	4.398	1.243	0.534	0.114		CSEO
beta2_pH	11	1	0.480	0.401	0.184	1.202	1.196	0.304	0.260		CSEO
beta0_pelagic	11		-0.201	-0.116	-1.202	0.143	1.680	0.312	-0.272		CSEO
beta1_pelagic	11		1.236	1.127	0.733	2.420	1.658	0.421	0.259		CSEO
beta3_pelagic	11		29.960	30.113	25.500	32.706	1.214	1.561	-0.098		CSEO
beta2_yellow	15		-8.445	-9.039	-19.868	-0.424	1.108	5.139	0.116		SSEI
beta0_pH	16	4	0.967	0.954	-0.150	2.117	1.111	0.572	0.023		SSEO
beta3_pH	16	4	41.413	41.504	40.262	42.460	2.191	0.612	-0.149		SSEO
beta0_pelagic	16		-0.515	-0.389	-1.433	0.046	1.638	0.422	-0.299		SSEO
beta1_pelagic	16		1.766	1.529	0.727	3.351	1.597	0.778	0.305		SSEO
beta2_pelagic	16		0.819	0.335	0.101	4.615	1.129	1.149	0.421		SSEO
mu_beta1_pH	2	3	3.157	0.122	0.000	15.672	1.379	4.511	0.673		NG
tau_beta1_pH	2	3	2.466	0.024	0.010	27.993	1.145	9.266	0.264		NG
mu_beta1_pH	3	3	0.004	0.000	0.000	0.047	1.205	0.011	0.353		PWSI
tau_beta1_pH	3	2	0.949	0.039	0.010	8.764	1.108	5.554	0.164		PWSI
tau_beta1_pH	3	4	5.886	3.111	0.301	31.577	1.104	10.010	0.277		PWSI
beta0_pelagic	3		0.213	0.252	-0.329	0.557	1.150	0.223	-0.175		PWSI
beta3_pelagic	3		28.080	28.498	23.552	32.866	1.126	2.688	-0.156		PWSI
beta0_pH	8	1	-0.568	-0.523	-1.377	-0.045	1.135	0.352	-0.128		afognak

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.01	615.00	211.00	1,011.00	339.00
Prop. of total	0.74	0.26	0.75	0.25
Prop. median less than 1	0.21	0.24	0.32	0.43
Prop. median less than 10	0.30	0.60	0.54	0.76
Prop. median less than 50	0.44	0.83	0.71	0.82
Median of the median values	91.35	5.98	6.49	1.82
Mean of the median values	758.32	138.35	197.42	147.41
Mean skew of posteriors	0.23	0.29	0.27	0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.1	1.00	12.00	56.00	48.00
Prop. of total	0.08	0.92	0.54	0.46
Prop. median less than 1	0.00	0.58	0.73	0.48
Prop. median less than 10	0.00	0.83	0.82	0.77
Prop. median less than 50	1.00	1.00	0.95	0.83
Median of the median values	16.12	0.47	0.01	1.58
Mean of the median values	16.12	9.64	35.55	20.34
Mean skew of posteriors	0.22	0.21	0.14	0.18

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable	Freq
beta0_pelagic	3
beta0_pH	3
beta1_pelagic	2
beta1_pH	1
beta2_pelagic	1
beta2_pH	1
beta3_pelagic	2
beta3_pH	1
H_ay	1
H_ayg	1
H_ayu	1
Hb_ayg	1
Hd_ayg	2
Ho_ay	1
Ho_ayg	1
Ho_ayu	1
Hp_ay	1
Hp_ayg	1
Hs_ay	2
Hs_ayg	2
Hs_ayu	5
Hy_ay	1
Hy_ayu	2
mu_beta1_pH	2
mu_wt	1
p_dsr	2
p_pelagic	28
p_slope	2
pH	12
R_ay	2
R_ayg	5
R_ayu	5
Rb_ay	1
Rb_ayg	4
Rb_ayu	3
Rd_ay	2
Rd_ayg	5
Rd_ayu	1
re_dsr	2
re_pelagic	14
re_pH	14
re_slope	2
Ro_ay	15
Ro_ayg	16
Ro_ayu	31
Rp_ay	1
Rp_ayg	4
Rp_ayu	3
Rs_ay	4
Rs_ayg	6
Rs_ayu	13
Ry_ay	12
Ry_ayg	6
Ry_ayu	25
tau_beta1_pH	1
wt	21

Parameter estimates:

Summary Table of Parameter Estimates

Parameter	mean	sd	Lower_CI	Median	Upper_CI
mu_bc_H[1]	-0.013	0.070	-0.142	-0.015	0.130
mu_bc_H[2]	0.069	0.053	-0.028	0.066	0.182
mu_bc_H[3]	-0.279	0.070	-0.414	-0.280	-0.140
mu_bc_H[4]	-0.927	0.186	-1.303	-0.919	-0.566
mu_bc_H[5]	1.307	1.137	0.083	1.052	4.226
mu_bc_H[6]	-2.431	0.318	-3.043	-2.436	-1.790
mu_bc_H[7]	-0.450	0.101	-0.653	-0.449	-0.259
mu_bc_H[8]	0.260	0.341	-0.313	0.234	1.038
mu_bc_H[9]	-0.247	0.127	-0.497	-0.250	0.001
mu_bc_H[10]	0.059	0.065	-0.063	0.057	0.190
mu_bc_H[11]	-0.126	0.037	-0.200	-0.125	-0.052
mu_bc_H[12]	-0.248	0.104	-0.470	-0.244	-0.047
mu_bc_H[13]	-0.160	0.073	-0.300	-0.161	-0.014
mu_bc_H[14]	-0.302	0.090	-0.486	-0.301	-0.133
mu_bc_H[15]	-0.337	0.047	-0.432	-0.336	-0.244
mu_bc_H[16]	-0.142	0.384	-0.774	-0.184	0.727
mu_bc_R[1]	1.053	0.123	0.811	1.051	1.287
mu_bc_R[2]	1.015	0.078	0.850	1.017	1.162
mu_bc_R[3]	0.237	0.093	0.048	0.239	0.416
mu_bc_R[4]	0.683	0.197	0.272	0.688	1.046
mu_bc_R[5]	1.681	0.380	0.928	1.673	2.429
mu_bc_R[6]	-1.457	0.337	-2.135	-1.463	-0.795
mu_bc_R[7]	0.514	0.167	0.177	0.516	0.830
mu_bc_R[8]	0.650	0.152	0.362	0.647	0.949
mu_bc_R[9]	0.395	0.170	0.054	0.401	0.723
mu_bc_R[10]	0.907	0.103	0.704	0.905	1.111
mu_bc_R[11]	1.110	0.076	0.964	1.110	1.266
mu_bc_R[12]	0.879	0.162	0.567	0.877	1.209
mu_bc_R[13]	1.036	0.093	0.849	1.037	1.214
mu_bc_R[14]	0.872	0.122	0.636	0.872	1.114
mu_bc_R[15]	0.802	0.097	0.611	0.803	0.988
mu_bc_R[16]	1.355	0.114	1.120	1.355	1.576
tau_pH[1]	2.418	0.183	2.070	2.419	2.783
tau_pH[2]	1.635	0.154	1.347	1.629	1.951
tau_pH[3]	1.206	0.112	1.000	1.202	1.437
tau_pH[4]	1.073	0.145	0.817	1.066	1.374
tau_pH[5]	0.445	0.079	0.305	0.438	0.616
beta0_pH[1,1]	0.581	0.178	0.199	0.588	0.899
beta0_pH[2,1]	1.412	0.164	1.084	1.417	1.732
beta0_pH[3,1]	1.401	0.207	0.922	1.419	1.751
beta0_pH[4,1]	1.573	0.250	0.975	1.613	1.966
beta0_pH[5,1]	-0.691	0.360	-1.519	-0.651	-0.113
beta0_pH[6,1]	-0.440	0.403	-1.493	-0.344	0.062
beta0_pH[7,1]	0.759	0.011	0.737	0.758	0.781
beta0_pH[8,1]	-0.568	0.352	-1.377	-0.523	-0.045
beta0_pH[9,1]	-0.286	0.275	-0.873	-0.259	0.197
beta0_pH[10,1]	0.304	0.197	-0.086	0.315	0.658
beta0_pH[11,1]	0.118	0.356	-0.610	0.144	0.696
beta0_pH[12,1]	0.614	0.196	0.209	0.623	0.966
beta0_pH[13,1]	-0.687	0.151	-1.008	-0.682	-0.411
beta0_pH[14,1]	0.009	0.168	-0.333	0.013	0.324
beta0_pH[15,1]	0.069	0.183	-0.330	0.081	0.385
beta0_pH[16,1]	-0.870	0.172	-1.172	-0.879	-0.533
beta0_pH[1,2]	2.563	0.085	2.403	2.559	2.736
beta0_pH[2,2]	3.079	0.022	3.036	3.079	3.123
beta0_pH[3,2]	3.625	0.040	3.549	3.625	3.702
beta0_pH[4,2]	2.936	0.019	2.899	2.936	2.974
beta0_pH[5,2]	4.536	1.438	2.783	4.268	7.985
beta0_pH[6,2]	3.322	0.158	3.027	3.319	3.646
beta0_pH[7,2]	2.034	0.049	1.938	2.034	2.132
beta0_pH[8,2]	2.822	0.064	2.698	2.822	2.950
beta0_pH[9,2]	3.952	0.206	3.574	3.948	4.376
beta0_pH[10,2]	3.993	0.130	3.748	3.989	4.258
beta0_pH[11,2]	1.885	0.009	1.868	1.885	1.902
beta0_pH[12,2]	1.793	0.068	1.661	1.793	1.931
beta0_pH[13,2]	2.497	0.027	2.446	2.498	2.548
beta0_pH[14,2]	1.803	0.024	1.756	1.802	1.851
beta0_pH[15,2]	2.545	0.017	2.511	2.545	2.578
beta0_pH[16,2]	2.765	0.025	2.715	2.766	2.814
beta0_pH[1,3]	1.085	0.427	0.108	1.172	1.650
beta0_pH[2,3]	2.125	0.017	2.092	2.125	2.159
beta0_pH[3,3]	2.691	0.024	2.644	2.691	2.737
beta0_pH[4,3]	2.980	0.035	2.913	2.980	3.049
beta0_pH[5,3]	3.312	2.164	0.944	2.777	9.082
beta0_pH[6,3]	1.365	0.373	0.718	1.340	2.203
beta0_pH[7,3]	0.860	0.050	0.763	0.860	0.961
beta0_pH[8,3]	0.274	0.340	-0.468	0.315	0.715
beta0_pH[9,3]	-2.076	1.109	-4.685	-1.956	-0.304
beta0_pH[10,3]	-0.333	0.775	-2.103	-0.270	0.956
beta0_pH[11,3]	1.537	0.538	0.478	1.532	2.598
beta0_pH[12,3]	0.950	0.798	-0.786	1.010	2.267
beta0_pH[13,3]	1.438	0.506	0.420	1.435	2.414
beta0_pH[14,3]	0.723	0.488	-0.343	0.761	1.567
beta0_pH[15,3]	0.995	0.484	0.037	0.991	1.986
beta0_pH[16,3]	1.182	0.643	-0.121	1.194	2.411
beta0_pH[11,4]	0.914	0.455	-0.005	0.898	1.900
beta0_pH[12,4]	0.840	0.717	-0.664	0.818	2.298
beta0_pH[13,4]	0.609	0.607	-0.888	0.701	1.579
beta0_pH[14,4]	0.359	0.442	-0.568	0.406	1.114
beta0_pH[15,4]	0.551	0.420	-0.581	0.588	1.262
beta0_pH[16,4]	0.967	0.572	-0.150	0.954	2.117
beta0_pH[11,5]	2.961	0.145	2.697	2.953	3.265
beta0_pH[12,5]	-0.917	0.504	-1.921	-0.912	0.037
beta0_pH[13,5]	4.862	0.420	4.096	4.852	5.753
beta0_pH[14,5]	3.881	0.359	3.220	3.867	4.648
beta0_pH[15,5]	2.285	0.123	2.048	2.282	2.542
beta0_pH[16,5]	0.259	0.067	0.163	0.247	0.423
beta1_pH[1,1]	2.510	0.463	1.765	2.463	3.596
beta1_pH[2,1]	1.606	0.308	1.038	1.599	2.178
beta1_pH[3,1]	1.804	0.392	1.171	1.760	2.761
beta1_pH[4,1]	1.977	0.566	1.220	1.850	3.644
beta1_pH[5,1]	2.177	0.449	1.435	2.128	3.211
beta1_pH[6,1]	1.897	0.605	1.060	1.781	3.533
beta1_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,1]	2.087	0.715	1.136	1.933	4.026
beta1_pH[9,1]	1.878	0.434	1.196	1.811	2.876
beta1_pH[10,1]	2.065	0.349	1.460	2.034	2.864
beta1_pH[11,1]	3.273	0.534	2.414	3.212	4.398
beta1_pH[12,1]	2.921	0.295	2.374	2.912	3.535
beta1_pH[13,1]	3.497	0.250	3.041	3.487	4.039
beta1_pH[14,1]	2.827	0.312	2.216	2.824	3.443
beta1_pH[15,1]	2.382	0.262	1.915	2.367	2.959
beta1_pH[16,1]	4.020	0.242	3.538	4.031	4.446
beta1_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[1,3]	1.951	1.236	0.112	1.787	5.180
beta1_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,3]	4.152	2.449	0.096	3.767	10.257
beta1_pH[9,3]	7.227	2.823	3.090	6.760	13.902
beta1_pH[10,3]	5.993	2.442	2.515	5.540	11.683
beta1_pH[11,3]	1.687	0.631	0.390	1.706	2.856
beta1_pH[12,3]	3.291	0.819	1.894	3.212	5.167
beta1_pH[13,3]	1.497	0.572	0.345	1.518	2.585
beta1_pH[14,3]	2.028	0.559	1.010	2.000	3.212
beta1_pH[15,3]	1.772	0.578	0.574	1.774	2.885
beta1_pH[16,3]	3.726	0.693	2.378	3.701	5.128
beta1_pH[11,4]	1.912	0.526	0.799	1.938	2.895
beta1_pH[12,4]	3.233	0.731	1.756	3.239	4.805
beta1_pH[13,4]	1.404	0.611	0.380	1.341	2.829
beta1_pH[14,4]	2.023	0.505	1.123	1.990	3.063
beta1_pH[15,4]	1.463	0.458	0.660	1.432	2.567
beta1_pH[16,4]	3.471	0.630	2.159	3.494	4.661
beta1_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,1]	1.322	0.878	0.328	1.030	3.379
beta2_pH[2,1]	1.628	0.962	0.269	1.459	3.431
beta2_pH[3,1]	1.346	0.833	0.236	1.179	2.939
beta2_pH[4,1]	1.073	0.787	0.167	0.815	2.892
beta2_pH[5,1]	2.604	2.717	0.328	1.730	10.613
beta2_pH[6,1]	1.765	2.807	0.142	0.651	10.714
beta2_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,1]	1.535	2.471	0.155	0.502	9.451
beta2_pH[9,1]	1.798	2.589	0.218	0.727	10.695
beta2_pH[10,1]	1.888	2.405	0.333	0.952	10.029
beta2_pH[11,1]	0.480	0.304	0.184	0.401	1.202
beta2_pH[12,1]	1.560	0.817	0.610	1.291	3.689
beta2_pH[13,1]	1.006	0.447	0.515	0.895	2.301
beta2_pH[14,1]	1.208	0.728	0.473	0.968	3.423
beta2_pH[15,1]	1.057	0.648	0.372	0.873	3.125
beta2_pH[16,1]	1.036	0.304	0.630	0.972	1.750
beta2_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,3]	7.070	6.214	0.194	5.673	20.133
beta2_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,3]	1.657	1.689	0.092	1.115	6.261
beta2_pH[9,3]	0.378	0.919	0.086	0.210	1.649
beta2_pH[10,3]	0.548	1.083	0.117	0.290	2.787
beta2_pH[11,3]	-7.247	5.112	-20.052	-5.954	-0.923
beta2_pH[12,3]	-7.731	4.941	-20.000	-6.600	-1.432
beta2_pH[13,3]	-7.375	4.960	-19.658	-6.083	-1.200
beta2_pH[14,3]	-7.355	4.986	-19.817	-6.005	-1.383
beta2_pH[15,3]	-7.464	5.044	-20.144	-6.116	-1.172
beta2_pH[16,3]	-7.200	4.903	-19.577	-5.776	-1.531
beta2_pH[11,4]	-6.061	4.427	-17.242	-4.863	-0.821
beta2_pH[12,4]	-6.306	4.524	-17.856	-5.097	-1.062
beta2_pH[13,4]	-6.110	4.477	-17.339	-4.874	-0.839
beta2_pH[14,4]	-6.205	4.392	-17.358	-5.069	-1.160
beta2_pH[15,4]	-6.221	4.435	-17.069	-5.056	-0.935
beta2_pH[16,4]	-6.243	4.215	-16.903	-5.031	-1.410
beta2_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta3_pH[1,1]	34.992	1.121	32.925	35.021	37.440
beta3_pH[2,1]	32.650	1.727	29.657	32.524	36.687
beta3_pH[3,1]	33.085	1.357	30.580	33.042	36.213
beta3_pH[4,1]	33.039	1.929	29.935	32.869	37.935
beta3_pH[5,1]	27.663	1.015	26.185	27.562	30.353
beta3_pH[6,1]	32.858	2.761	25.813	32.976	38.241
beta3_pH[7,1]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,1]	34.426	3.442	29.165	33.807	41.099
beta3_pH[9,1]	31.768	2.295	27.779	31.685	36.588
beta3_pH[10,1]	32.458	1.344	30.011	32.393	35.534
beta3_pH[11,1]	31.414	1.193	29.373	31.375	33.872
beta3_pH[12,1]	30.779	0.469	29.814	30.801	31.636
beta3_pH[13,1]	31.376	0.539	30.348	31.362	32.437
beta3_pH[14,1]	32.413	0.747	31.147	32.350	33.958
beta3_pH[15,1]	30.988	0.724	29.487	31.020	32.359
beta3_pH[16,1]	30.107	0.352	29.429	30.095	30.787
beta3_pH[1,2]	34.415	5.526	25.530	34.345	43.539
beta3_pH[2,2]	34.706	5.465	25.530	34.684	43.555
beta3_pH[3,2]	30.449	7.270	18.748	30.333	42.423
beta3_pH[4,2]	30.280	7.267	18.466	30.159	42.337
beta3_pH[5,2]	30.712	7.269	18.553	30.914	42.403
beta3_pH[6,2]	30.107	7.209	18.513	30.083	42.261
beta3_pH[7,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[9,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[10,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[11,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[12,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[13,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[14,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[15,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[16,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[1,3]	38.432	3.060	31.413	39.574	42.823
beta3_pH[2,3]	34.593	5.549	25.393	34.710	43.633
beta3_pH[3,3]	30.731	7.180	18.690	30.791	42.467
beta3_pH[4,3]	30.398	7.147	18.568	30.296	42.194
beta3_pH[5,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[6,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[7,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,3]	41.746	1.548	36.915	41.852	43.826
beta3_pH[9,3]	36.210	2.014	34.522	35.402	42.209
beta3_pH[10,3]	37.641	2.409	34.640	36.965	43.298
beta3_pH[11,3]	40.383	1.194	37.994	40.366	42.561
beta3_pH[12,3]	42.005	0.757	40.280	42.261	42.858
beta3_pH[13,3]	41.105	0.884	39.230	41.160	42.670
beta3_pH[14,3]	40.385	0.728	38.956	40.405	41.725
beta3_pH[15,3]	40.476	0.842	38.290	40.566	41.890
beta3_pH[16,3]	41.826	0.385	41.023	41.893	42.559
beta3_pH[11,4]	40.337	1.079	38.175	40.281	42.539
beta3_pH[12,4]	41.855	0.795	40.239	42.128	42.892
beta3_pH[13,4]	41.789	1.251	39.206	41.823	43.735
beta3_pH[14,4]	40.445	0.699	39.018	40.477	41.749
beta3_pH[15,4]	41.303	0.985	39.469	41.213	43.382
beta3_pH[16,4]	41.413	0.612	40.262	41.504	42.460
beta3_pH[11,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[12,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[13,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[14,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[15,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[16,5]	25.000	0.000	25.000	25.000	25.000
beta4_pH[1,1]	0.264	0.021	0.224	0.264	0.306
beta4_pH[2,1]	0.211	0.005	0.200	0.211	0.221
beta4_pH[3,1]	0.173	0.007	0.158	0.173	0.186
beta4_pH[4,1]	0.293	0.011	0.271	0.294	0.316
beta4_pH[5,1]	0.423	0.264	0.088	0.388	0.907
beta4_pH[6,1]	0.786	0.327	0.330	0.747	1.505
beta4_pH[7,1]	0.846	0.332	0.353	0.794	1.670
beta4_pH[8,1]	0.160	0.071	0.045	0.152	0.320
beta4_pH[9,1]	0.210	0.085	0.062	0.206	0.388
beta4_pH[10,1]	0.001	0.001	0.000	0.001	0.003
beta4_pH[11,1]	0.447	0.012	0.424	0.447	0.470
beta4_pH[12,1]	0.641	0.022	0.599	0.640	0.687
beta4_pH[13,1]	0.630	0.014	0.603	0.630	0.658
beta4_pH[14,1]	0.669	0.030	0.612	0.669	0.732
beta4_pH[15,1]	0.798	0.021	0.758	0.797	0.840
beta4_pH[16,1]	1.221	0.041	1.144	1.220	1.303
beta4_pH[1,2]	0.306	0.095	0.132	0.300	0.502
beta4_pH[2,2]	0.513	0.022	0.470	0.514	0.558
beta4_pH[3,2]	0.446	0.012	0.423	0.446	0.469
beta4_pH[4,2]	0.654	0.038	0.580	0.654	0.730
beta4_pH[5,2]	0.727	0.761	0.119	0.544	2.416
beta4_pH[6,2]	0.935	0.893	0.198	0.748	2.812
beta4_pH[7,2]	0.893	0.409	0.348	0.817	1.898
beta4_pH[8,2]	0.852	0.199	0.490	0.841	1.256
beta4_pH[9,2]	1.197	1.296	0.307	0.892	3.978
beta4_pH[10,2]	0.456	0.052	0.357	0.456	0.561
beta4_pH[11,2]	1.271	0.051	1.177	1.270	1.373
beta4_pH[12,2]	0.578	0.109	0.380	0.573	0.813
beta4_pH[13,2]	0.945	0.045	0.859	0.944	1.033
beta4_pH[14,2]	0.930	0.061	0.814	0.928	1.053
beta4_pH[15,2]	1.010	0.025	0.963	1.010	1.061
beta4_pH[16,2]	1.409	0.075	1.265	1.407	1.563
beta4_pH[1,3]	0.283	0.086	0.125	0.278	0.467
beta4_pH[2,3]	0.051	0.017	0.020	0.051	0.088
beta4_pH[3,3]	0.052	0.010	0.034	0.052	0.071
beta4_pH[4,3]	0.173	0.032	0.112	0.172	0.235
beta4_pH[5,3]	0.747	0.782	0.125	0.558	2.540
beta4_pH[6,3]	0.879	0.782	0.197	0.724	2.558
beta4_pH[7,3]	0.896	0.475	0.327	0.813	1.897
beta4_pH[8,3]	0.418	0.198	0.108	0.391	0.869
beta4_pH[9,3]	0.932	0.855	0.161	0.724	3.010
beta4_pH[10,3]	0.061	0.041	0.003	0.055	0.157
beta4_pH[11,3]	0.364	0.016	0.334	0.364	0.394
beta4_pH[12,3]	0.439	0.041	0.363	0.438	0.518
beta4_pH[13,3]	0.840	0.037	0.766	0.839	0.913
beta4_pH[14,3]	0.823	0.066	0.702	0.819	0.963
beta4_pH[15,3]	0.636	0.016	0.605	0.635	0.670
beta4_pH[16,3]	0.634	0.034	0.570	0.634	0.705
beta4_pH[1,4]	0.641	0.352	0.212	0.569	1.483
beta4_pH[2,4]	0.825	0.853	0.113	0.598	2.836
beta4_pH[3,4]	0.868	1.017	0.112	0.584	3.176
beta4_pH[4,4]	0.740	0.493	0.214	0.630	1.908
beta4_pH[5,4]	0.740	0.778	0.116	0.556	2.489
beta4_pH[6,4]	0.916	0.773	0.219	0.737	2.598
beta4_pH[7,4]	0.900	0.430	0.363	0.821	1.925
beta4_pH[8,4]	0.749	0.656	0.161	0.592	2.180
beta4_pH[9,4]	0.779	0.825	0.123	0.574	2.640
beta4_pH[10,4]	1.448	3.019	0.007	0.389	11.109
beta4_pH[11,4]	0.444	0.017	0.411	0.444	0.478
beta4_pH[12,4]	0.557	0.052	0.459	0.556	0.666
beta4_pH[13,4]	0.925	0.037	0.853	0.925	0.996
beta4_pH[14,4]	0.918	0.074	0.777	0.915	1.074
beta4_pH[15,4]	0.713	0.020	0.675	0.713	0.753
beta4_pH[16,4]	0.688	0.041	0.610	0.687	0.770
beta4_pH[1,5]	0.625	0.347	0.218	0.553	1.481
beta4_pH[2,5]	0.804	0.828	0.108	0.582	3.029
beta4_pH[3,5]	0.858	0.943	0.115	0.587	3.335
beta4_pH[4,5]	0.732	0.447	0.218	0.623	1.915
beta4_pH[5,5]	0.768	0.800	0.117	0.558	2.689
beta4_pH[6,5]	0.881	0.777	0.195	0.728	2.546
beta4_pH[7,5]	0.900	0.453	0.344	0.823	1.901
beta4_pH[8,5]	0.743	0.617	0.151	0.603	2.116
beta4_pH[9,5]	0.777	0.844	0.124	0.563	2.694
beta4_pH[10,5]	1.479	3.068	0.006	0.392	10.771
beta4_pH[11,5]	0.603	0.051	0.507	0.602	0.704
beta4_pH[12,5]	0.586	0.166	0.306	0.570	0.965
beta4_pH[13,5]	0.778	0.061	0.664	0.776	0.901
beta4_pH[14,5]	0.840	0.117	0.621	0.837	1.080
beta4_pH[15,5]	1.143	0.073	1.005	1.141	1.293
beta4_pH[16,5]	18.173	4.083	10.100	18.353	24.624
beta0_pelagic[1]	2.269	0.036	2.202	2.268	2.342
beta0_pelagic[2]	1.538	0.007	1.524	1.538	1.551
beta0_pelagic[3]	0.213	0.223	-0.329	0.252	0.557
beta0_pelagic[4]	0.186	0.235	-0.413	0.210	0.597
beta0_pelagic[5]	1.450	0.154	1.153	1.446	1.758
beta0_pelagic[6]	1.355	0.060	1.242	1.354	1.476
beta0_pelagic[7]	1.484	0.039	1.411	1.483	1.562
beta0_pelagic[8]	1.674	0.043	1.592	1.673	1.763
beta0_pelagic[9]	2.795	0.066	2.671	2.794	2.925
beta0_pelagic[10]	2.602	0.032	2.540	2.601	2.667
beta0_pelagic[11]	-0.201	0.312	-1.202	-0.116	0.143
beta0_pelagic[12]	1.800	0.036	1.735	1.799	1.876
beta0_pelagic[13]	-0.691	0.372	-1.408	-0.654	-0.011
beta0_pelagic[14]	-0.277	0.179	-0.691	-0.258	0.019
beta0_pelagic[15]	-0.221	0.040	-0.321	-0.204	-0.182
beta0_pelagic[16]	-0.515	0.422	-1.433	-0.389	0.046
beta1_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[3]	0.740	0.245	0.329	0.718	1.369
beta1_pelagic[4]	0.874	0.275	0.417	0.849	1.551
beta1_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[11]	1.236	0.421	0.733	1.127	2.420
beta1_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[13]	2.392	0.579	1.294	2.376	3.420
beta1_pelagic[14]	1.805	0.479	1.110	1.714	2.998
beta1_pelagic[15]	0.160	0.217	0.001	0.044	0.730
beta1_pelagic[16]	1.766	0.778	0.727	1.529	3.351
beta2_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[3]	6.770	5.975	0.207	5.288	21.086
beta2_pelagic[4]	7.449	6.155	0.223	5.999	21.932
beta2_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[11]	1.683	1.510	0.167	1.219	5.796
beta2_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[13]	0.234	0.241	0.104	0.191	0.556
beta2_pelagic[14]	0.895	0.981	0.202	0.572	3.760
beta2_pelagic[15]	2.684	1.858	0.188	2.363	7.109
beta2_pelagic[16]	0.819	1.149	0.101	0.335	4.615
beta3_pelagic[1]	29.619	3.729	23.325	29.680	35.715
beta3_pelagic[2]	29.540	3.764	23.348	29.602	35.692
beta3_pelagic[3]	28.080	2.688	23.552	28.498	32.866
beta3_pelagic[4]	26.086	1.344	23.615	26.236	28.777
beta3_pelagic[5]	29.641	6.678	18.705	29.641	40.635
beta3_pelagic[6]	29.743	6.746	18.578	29.734	40.511
beta3_pelagic[7]	29.656	6.564	18.627	29.969	40.422
beta3_pelagic[8]	29.424	6.665	18.621	29.316	40.339
beta3_pelagic[9]	29.620	6.713	18.571	29.721	40.461
beta3_pelagic[10]	29.361	6.661	18.505	29.362	40.336
beta3_pelagic[11]	29.960	1.561	25.500	30.113	32.706
beta3_pelagic[12]	32.374	5.509	23.388	32.305	41.634
beta3_pelagic[13]	32.192	3.558	26.293	31.679	39.945
beta3_pelagic[14]	34.355	1.733	31.632	34.096	38.591
beta3_pelagic[15]	35.039	5.753	23.682	36.353	41.805
beta3_pelagic[16]	31.717	3.262	25.431	31.296	39.447
beta0_yellow[1]	-0.483	0.079	-0.643	-0.483	-0.330
beta0_yellow[2]	0.543	0.080	0.377	0.547	0.683
beta0_yellow[3]	-0.202	0.130	-0.459	-0.205	0.018
beta0_yellow[4]	0.926	0.203	0.456	0.984	1.184
beta0_yellow[5]	0.256	0.330	-0.448	0.278	0.839
beta0_yellow[6]	0.818	0.105	0.608	0.820	1.025
beta0_yellow[7]	1.000	0.094	0.816	1.001	1.183
beta0_yellow[8]	1.022	0.094	0.840	1.021	1.209
beta0_yellow[9]	0.458	0.101	0.264	0.459	0.657
beta0_yellow[10]	0.538	0.071	0.399	0.537	0.681
beta0_yellow[11]	-0.195	0.121	-0.454	-0.185	0.017
beta0_yellow[12]	-1.087	0.094	-1.269	-1.087	-0.904
beta0_yellow[13]	-1.364	0.208	-1.750	-1.345	-1.066
beta0_yellow[14]	-0.257	0.174	-0.600	-0.235	-0.015
beta0_yellow[15]	-1.251	0.133	-1.563	-1.236	-1.028
beta0_yellow[16]	-0.254	0.094	-0.441	-0.252	-0.069
beta1_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[2]	1.012	0.238	0.571	1.000	1.489
beta1_yellow[3]	0.451	0.240	0.020	0.463	0.896
beta1_yellow[4]	1.249	0.492	0.635	1.155	2.746
beta1_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[11]	0.313	0.387	0.001	0.128	1.252
beta1_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[13]	1.136	1.266	0.023	0.775	5.578
beta1_yellow[14]	0.753	1.089	0.013	0.431	4.552
beta1_yellow[15]	0.185	0.228	0.011	0.092	0.769
beta1_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[2]	-7.685	5.842	-21.611	-6.326	-0.465
beta2_yellow[3]	-7.670	5.870	-21.780	-6.404	-0.411
beta2_yellow[4]	-5.950	5.861	-20.715	-4.179	-0.154
beta2_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[11]	-7.739	5.986	-22.377	-6.294	-0.423
beta2_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[13]	-7.505	6.364	-21.262	-5.932	-0.320
beta2_yellow[14]	-7.317	5.823	-21.486	-5.905	-0.372
beta2_yellow[15]	-8.445	5.139	-19.868	-9.039	-0.424
beta2_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta3_yellow[1]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[2]	28.854	1.124	27.258	28.680	31.721
beta3_yellow[3]	32.353	3.492	22.932	32.601	39.299
beta3_yellow[4]	29.941	3.210	24.615	28.704	35.746
beta3_yellow[5]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[6]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[7]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[8]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[9]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[10]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[11]	34.366	3.472	29.474	33.459	42.278
beta3_yellow[12]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[13]	33.651	2.614	29.298	33.812	39.463
beta3_yellow[14]	33.785	3.577	29.184	33.248	43.084
beta3_yellow[15]	35.775	3.620	30.237	35.625	43.360
beta3_yellow[16]	25.000	0.000	25.000	25.000	25.000
beta0_black[1]	0.012	0.075	-0.132	0.012	0.161
beta0_black[2]	1.855	0.032	1.792	1.855	1.920
beta0_black[3]	1.239	0.037	1.170	1.239	1.313
beta0_black[4]	2.425	0.047	2.336	2.424	2.522
beta0_black[5]	5.054	2.282	1.870	4.616	10.699
beta0_black[6]	5.165	2.123	2.617	4.593	10.724
beta0_black[7]	4.102	2.109	1.618	3.541	9.922
beta0_black[8]	0.878	0.216	0.479	0.874	1.332
beta0_black[9]	2.815	0.274	2.280	2.810	3.361
beta0_black[10]	1.288	0.051	1.185	1.288	1.387
beta0_black[11]	2.978	0.051	2.879	2.978	3.078
beta0_black[12]	4.918	0.190	4.551	4.917	5.289
beta0_black[13]	0.442	0.327	0.019	0.377	1.214
beta0_black[14]	2.851	0.107	2.644	2.850	3.062
beta0_black[15]	0.723	0.051	0.622	0.723	0.823
beta0_black[16]	3.772	0.073	3.633	3.770	3.916
beta2_black[1]	0.000	0.000	0.000	0.000	0.000
beta2_black[2]	0.000	0.000	0.000	0.000	0.000
beta2_black[3]	0.000	0.000	0.000	0.000	0.000
beta2_black[4]	0.000	0.000	0.000	0.000	0.000
beta2_black[5]	0.000	0.000	0.000	0.000	0.000
beta2_black[6]	0.000	0.000	0.000	0.000	0.000
beta2_black[7]	0.000	0.000	0.000	0.000	0.000
beta2_black[8]	0.000	0.000	0.000	0.000	0.000
beta2_black[9]	0.000	0.000	0.000	0.000	0.000
beta2_black[10]	0.000	0.000	0.000	0.000	0.000
beta2_black[11]	0.000	0.000	0.000	0.000	0.000
beta2_black[12]	0.000	0.000	0.000	0.000	0.000
beta2_black[13]	-2.534	1.738	-6.801	-2.149	-0.295
beta2_black[14]	0.000	0.000	0.000	0.000	0.000
beta2_black[15]	0.000	0.000	0.000	0.000	0.000
beta2_black[16]	0.000	0.000	0.000	0.000	0.000
beta3_black[1]	25.000	0.000	25.000	25.000	25.000
beta3_black[2]	25.000	0.000	25.000	25.000	25.000
beta3_black[3]	25.000	0.000	25.000	25.000	25.000
beta3_black[4]	25.000	0.000	25.000	25.000	25.000
beta3_black[5]	25.000	0.000	25.000	25.000	25.000
beta3_black[6]	25.000	0.000	25.000	25.000	25.000
beta3_black[7]	25.000	0.000	25.000	25.000	25.000
beta3_black[8]	25.000	0.000	25.000	25.000	25.000
beta3_black[9]	25.000	0.000	25.000	25.000	25.000
beta3_black[10]	25.000	0.000	25.000	25.000	25.000
beta3_black[11]	25.000	0.000	25.000	25.000	25.000
beta3_black[12]	25.000	0.000	25.000	25.000	25.000
beta3_black[13]	37.448	3.200	27.433	38.201	40.874
beta3_black[14]	25.000	0.000	25.000	25.000	25.000
beta3_black[15]	25.000	0.000	25.000	25.000	25.000
beta3_black[16]	25.000	0.000	25.000	25.000	25.000
beta4_black[1]	-0.326	0.091	-0.505	-0.326	-0.145
beta4_black[2]	0.258	0.053	0.155	0.259	0.362
beta4_black[3]	-0.719	0.064	-0.842	-0.721	-0.595
beta4_black[4]	0.628	0.132	0.380	0.628	0.886
beta4_black[5]	0.767	1.527	-1.183	0.475	4.304
beta4_black[6]	0.826	1.307	-0.876	0.554	4.030
beta4_black[7]	0.738	1.320	-0.909	0.492	3.696
beta4_black[8]	0.080	0.284	-0.520	0.089	0.623
beta4_black[9]	1.019	0.824	-0.109	0.858	3.118
beta4_black[10]	0.074	0.067	-0.055	0.072	0.210
beta4_black[11]	-0.608	0.067	-0.737	-0.608	-0.473
beta4_black[12]	1.239	0.294	0.683	1.228	1.823
beta4_black[13]	-1.004	0.076	-1.153	-1.005	-0.855
beta4_black[14]	0.221	0.180	-0.126	0.222	0.590
beta4_black[15]	0.002	0.046	-0.085	0.003	0.093
beta4_black[16]	-0.415	0.116	-0.647	-0.414	-0.186
beta0_dsr[11]	-0.328	0.009	-0.347	-0.328	-0.311
beta0_dsr[12]	1.408	0.045	1.320	1.407	1.498
beta0_dsr[13]	0.155	0.022	0.113	0.155	0.196
beta0_dsr[14]	-0.127	0.017	-0.160	-0.127	-0.095
beta0_dsr[15]	0.206	0.014	0.180	0.206	0.232
beta0_dsr[16]	0.160	0.016	0.129	0.160	0.191
beta1_dsr[11]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[12]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[13]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[14]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[15]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[16]	0.000	0.000	0.000	0.000	0.000
beta4_dsr[11]	0.545	0.035	0.473	0.545	0.614
beta4_dsr[12]	0.798	0.230	0.356	0.791	1.264
beta4_dsr[13]	-0.077	0.040	-0.154	-0.077	0.002
beta4_dsr[14]	-0.208	0.071	-0.346	-0.208	-0.071
beta4_dsr[15]	0.216	0.021	0.175	0.216	0.258
beta4_dsr[16]	-0.050	0.045	-0.138	-0.050	0.037
beta0_slope[11]	-2.209	0.021	-2.250	-2.209	-2.168
beta0_slope[12]	-4.984	0.188	-5.371	-4.976	-4.636
beta0_slope[13]	-1.501	0.035	-1.570	-1.502	-1.435
beta0_slope[14]	-2.712	0.051	-2.812	-2.710	-2.614
beta0_slope[15]	-1.592	0.022	-1.635	-1.592	-1.550
beta0_slope[16]	-2.657	0.040	-2.733	-2.658	-2.577
beta1_slope[11]	0.000	0.000	0.000	0.000	0.000
beta1_slope[12]	0.000	0.000	0.000	0.000	0.000
beta1_slope[13]	0.000	0.000	0.000	0.000	0.000
beta1_slope[14]	0.000	0.000	0.000	0.000	0.000
beta1_slope[15]	0.000	0.000	0.000	0.000	0.000
beta1_slope[16]	0.000	0.000	0.000	0.000	0.000
beta4_slope[11]	-0.334	0.074	-0.488	-0.333	-0.190
beta4_slope[12]	-4.652	0.791	-6.323	-4.620	-3.209
beta4_slope[13]	0.327	0.054	0.225	0.326	0.436
beta4_slope[14]	-0.157	0.153	-0.466	-0.155	0.133
beta4_slope[15]	-0.360	0.033	-0.424	-0.361	-0.295
beta4_slope[16]	-0.134	0.103	-0.346	-0.132	0.062
beta5_slope[11]	3.067	0.108	2.865	3.068	3.289
beta5_slope[12]	2.009	0.386	1.285	1.997	2.819
beta5_slope[13]	0.880	0.079	0.726	0.881	1.029
beta5_slope[14]	3.746	0.346	3.130	3.734	4.467
beta5_slope[15]	1.310	0.072	1.177	1.310	1.452
beta5_slope[16]	2.724	0.236	2.273	2.714	3.212
tau_comp[1]	3.790	0.298	3.239	3.779	4.430
tau_comp[2]	2.866	0.249	2.416	2.851	3.376
tau_comp[3]	0.870	0.074	0.737	0.868	1.029
tau_comp[4]	1.783	0.276	1.307	1.760	2.375
tau_comp[5]	1.280	0.187	0.948	1.264	1.677
sigma_H[1]	0.206	0.049	0.118	0.203	0.310
sigma_H[2]	0.174	0.028	0.127	0.172	0.238
sigma_H[3]	0.188	0.041	0.114	0.184	0.274
sigma_H[4]	0.462	0.081	0.330	0.453	0.644
sigma_H[5]	0.961	0.202	0.593	0.951	1.396
sigma_H[6]	0.552	0.193	0.189	0.540	0.953
sigma_H[7]	0.310	0.060	0.215	0.302	0.447
sigma_H[8]	0.425	0.088	0.298	0.413	0.618
sigma_H[9]	0.511	0.117	0.325	0.498	0.781
sigma_H[10]	0.211	0.039	0.145	0.208	0.300
sigma_H[11]	0.266	0.045	0.190	0.262	0.363
sigma_H[12]	0.386	0.149	0.196	0.340	0.726
sigma_H[13]	0.219	0.037	0.156	0.216	0.302
sigma_H[14]	0.498	0.091	0.343	0.492	0.694
sigma_H[15]	0.231	0.038	0.166	0.227	0.316
sigma_H[16]	0.214	0.040	0.144	0.210	0.306
lambda_H[1]	2.152	3.477	0.072	0.929	11.507
lambda_H[2]	5.205	4.974	0.508	3.667	19.039
lambda_H[3]	5.525	7.803	0.283	2.980	26.653
lambda_H[4]	0.006	0.004	0.001	0.005	0.017
lambda_H[5]	2.386	5.273	0.025	0.641	17.829
lambda_H[6]	6.693	14.088	0.009	0.451	47.082
lambda_H[7]	0.012	0.008	0.002	0.010	0.034
lambda_H[8]	7.711	10.210	0.087	4.237	36.462
lambda_H[9]	0.015	0.010	0.003	0.013	0.042
lambda_H[10]	0.244	0.662	0.025	0.152	0.848
lambda_H[11]	0.233	0.330	0.010	0.114	1.043
lambda_H[12]	3.818	5.326	0.171	2.108	17.935
lambda_H[13]	3.132	2.854	0.264	2.282	11.047
lambda_H[14]	3.280	3.963	0.243	2.021	13.658
lambda_H[15]	0.018	0.021	0.003	0.014	0.055
lambda_H[16]	0.486	0.725	0.034	0.252	2.361
mu_lambda_H[1]	4.108	1.905	1.093	3.904	8.301
mu_lambda_H[2]	3.708	1.917	0.514	3.547	7.620
mu_lambda_H[3]	3.439	1.919	0.694	3.137	7.743
sigma_lambda_H[1]	8.121	4.317	1.731	7.363	17.876
sigma_lambda_H[2]	8.211	4.666	0.850	7.590	18.204
sigma_lambda_H[3]	6.358	4.225	0.906	5.445	17.022
beta_H[1,1]	6.483	1.314	3.276	6.687	8.438
beta_H[2,1]	9.752	0.594	8.407	9.781	10.824
beta_H[3,1]	7.833	0.809	5.948	7.922	9.095
beta_H[4,1]	10.806	7.626	-4.443	10.611	25.915
beta_H[5,1]	-0.252	2.640	-5.707	-0.102	4.383
beta_H[6,1]	3.221	4.371	-7.612	4.841	8.203
beta_H[7,1]	1.221	5.679	-10.433	1.535	11.484
beta_H[8,1]	1.289	3.920	-2.743	1.214	3.445
beta_H[9,1]	13.615	5.671	2.437	13.466	25.136
beta_H[10,1]	6.820	1.815	2.887	6.860	10.352
beta_H[11,1]	4.980	3.603	-3.051	5.687	10.032
beta_H[12,1]	2.601	1.090	0.720	2.520	4.993
beta_H[13,1]	9.078	0.890	7.162	9.144	10.554
beta_H[14,1]	2.129	1.037	0.009	2.144	4.120
beta_H[15,1]	-7.168	3.270	-13.278	-7.296	-0.243
beta_H[16,1]	3.982	3.017	-0.782	3.497	11.041
beta_H[1,2]	7.776	0.273	7.232	7.780	8.307
beta_H[2,2]	9.884	0.161	9.560	9.884	10.207
beta_H[3,2]	8.798	0.204	8.409	8.796	9.206
beta_H[4,2]	2.983	1.460	0.183	2.943	5.940
beta_H[5,2]	1.775	0.976	-0.113	1.782	3.685
beta_H[6,2]	6.100	1.083	3.613	6.261	7.811
beta_H[7,2]	2.412	1.051	0.474	2.362	4.527
beta_H[8,2]	3.005	1.087	1.371	3.140	4.190
beta_H[9,2]	3.152	1.102	0.988	3.128	5.371
beta_H[10,2]	8.027	0.366	7.280	8.036	8.722
beta_H[11,2]	9.809	0.642	8.866	9.688	11.238
beta_H[12,2]	3.942	0.370	3.264	3.927	4.696
beta_H[13,2]	9.148	0.245	8.686	9.143	9.652
beta_H[14,2]	4.017	0.340	3.369	4.013	4.709
beta_H[15,2]	11.577	0.594	10.312	11.604	12.688
beta_H[16,2]	4.464	0.827	2.955	4.448	6.113
beta_H[1,3]	8.448	0.265	7.955	8.439	8.990
beta_H[2,3]	9.874	0.145	9.599	9.870	10.165
beta_H[3,3]	9.477	0.161	9.168	9.472	9.796
beta_H[4,3]	-2.128	0.875	-3.870	-2.121	-0.477
beta_H[5,3]	3.920	0.602	2.722	3.928	5.068
beta_H[6,3]	8.383	1.190	6.730	8.045	10.923
beta_H[7,3]	-2.584	0.619	-3.779	-2.582	-1.396
beta_H[8,3]	5.295	0.529	4.689	5.223	6.298
beta_H[9,3]	-2.666	0.694	-4.017	-2.662	-1.333
beta_H[10,3]	8.566	0.282	8.019	8.565	9.135
beta_H[11,3]	8.509	0.293	7.891	8.534	9.019
beta_H[12,3]	5.261	0.317	4.511	5.302	5.765
beta_H[13,3]	8.855	0.174	8.511	8.859	9.194
beta_H[14,3]	5.776	0.269	5.174	5.801	6.256
beta_H[15,3]	10.268	0.286	9.694	10.263	10.831
beta_H[16,3]	5.900	0.582	4.732	5.916	6.931
beta_H[1,4]	8.093	0.184	7.702	8.106	8.428
beta_H[2,4]	9.857	0.128	9.589	9.864	10.094
beta_H[3,4]	9.972	0.152	9.645	9.981	10.250
beta_H[4,4]	11.803	0.467	10.868	11.798	12.730
beta_H[5,4]	5.610	0.750	4.438	5.518	7.333
beta_H[6,4]	7.272	0.926	5.289	7.550	8.593
beta_H[7,4]	8.539	0.341	7.843	8.541	9.228
beta_H[8,4]	6.709	0.258	6.221	6.723	7.126
beta_H[9,4]	7.420	0.442	6.541	7.409	8.304
beta_H[10,4]	7.526	0.225	7.111	7.516	7.996
beta_H[11,4]	9.509	0.195	9.116	9.509	9.884
beta_H[12,4]	7.249	0.201	6.870	7.249	7.668
beta_H[13,4]	9.065	0.142	8.778	9.064	9.334
beta_H[14,4]	7.833	0.214	7.426	7.834	8.249
beta_H[15,4]	9.478	0.221	9.048	9.478	9.898
beta_H[16,4]	9.430	0.240	8.979	9.420	9.899
beta_H[1,5]	8.888	0.156	8.562	8.891	9.183
beta_H[2,5]	10.676	0.100	10.493	10.674	10.889
beta_H[3,5]	10.763	0.162	10.466	10.756	11.089
beta_H[4,5]	8.346	0.487	7.407	8.323	9.372
beta_H[5,5]	5.175	0.605	3.750	5.252	6.129
beta_H[6,5]	8.995	0.635	8.059	8.860	10.442
beta_H[7,5]	6.577	0.338	5.922	6.563	7.253
beta_H[8,5]	8.243	0.215	7.895	8.227	8.687
beta_H[9,5]	8.195	0.455	7.294	8.194	9.094
beta_H[10,5]	10.036	0.231	9.579	10.038	10.492
beta_H[11,5]	11.585	0.219	11.165	11.581	12.035
beta_H[12,5]	8.612	0.191	8.244	8.612	8.989
beta_H[13,5]	10.140	0.131	9.895	10.135	10.410
beta_H[14,5]	9.241	0.214	8.835	9.234	9.685
beta_H[15,5]	11.309	0.236	10.853	11.303	11.774
beta_H[16,5]	9.824	0.185	9.445	9.834	10.175
beta_H[1,6]	10.131	0.249	9.736	10.096	10.759
beta_H[2,6]	11.471	0.117	11.243	11.470	11.704
beta_H[3,6]	10.619	0.154	10.288	10.628	10.899
beta_H[4,6]	12.916	0.889	11.120	12.940	14.628
beta_H[5,6]	5.668	0.642	4.528	5.637	6.984
beta_H[6,6]	8.768	0.580	7.435	8.808	9.809
beta_H[7,6]	10.535	0.581	9.392	10.555	11.662
beta_H[8,6]	9.522	0.281	9.004	9.540	9.967
beta_H[9,6]	8.282	0.769	6.833	8.262	9.839
beta_H[10,6]	9.142	0.353	8.383	9.163	9.762
beta_H[11,6]	10.705	0.364	9.939	10.733	11.333
beta_H[12,6]	9.382	0.238	8.934	9.384	9.885
beta_H[13,6]	11.100	0.152	10.802	11.095	11.414
beta_H[14,6]	9.912	0.282	9.341	9.910	10.469
beta_H[15,6]	10.553	0.422	9.704	10.548	11.385
beta_H[16,6]	10.596	0.257	10.054	10.612	11.093
beta_H[1,7]	10.123	1.299	6.729	10.433	11.927
beta_H[2,7]	12.104	0.516	11.013	12.137	13.035
beta_H[3,7]	10.303	0.656	8.853	10.352	11.433
beta_H[4,7]	1.353	4.612	-7.366	1.332	10.884
beta_H[5,7]	5.853	2.093	1.874	5.908	9.967
beta_H[6,7]	7.352	3.181	-0.782	8.415	11.269
beta_H[7,7]	7.679	2.930	1.917	7.656	13.528
beta_H[8,7]	10.899	1.049	9.193	10.876	12.662
beta_H[9,7]	4.934	3.987	-3.242	4.948	12.484
beta_H[10,7]	10.020	1.656	7.142	9.885	13.712
beta_H[11,7]	11.303	1.782	8.282	11.104	15.202
beta_H[12,7]	9.794	0.996	7.770	9.879	11.454
beta_H[13,7]	11.716	0.712	10.014	11.782	12.913
beta_H[14,7]	10.591	0.900	8.752	10.635	12.325
beta_H[15,7]	13.981	2.177	9.832	13.959	18.371
beta_H[16,7]	12.833	1.376	10.501	12.702	15.845
beta0_H[1]	8.730	18.323	-28.593	8.560	47.155
beta0_H[2]	10.596	7.938	-5.558	10.603	26.246
beta0_H[3]	9.740	9.406	-9.186	9.904	29.320
beta0_H[4]	10.280	180.351	-347.308	8.680	374.330
beta0_H[5]	3.927	30.610	-56.011	3.928	66.200
beta0_H[6]	6.711	52.085	-108.424	7.618	117.609
beta0_H[7]	1.569	137.360	-297.040	7.808	269.865
beta0_H[8]	6.262	33.813	-20.610	6.455	32.601
beta0_H[9]	5.379	122.721	-234.777	3.659	252.161
beta0_H[10]	8.171	37.253	-71.961	8.591	78.070
beta0_H[11]	10.300	53.177	-91.563	9.395	125.753
beta0_H[12]	6.935	12.537	-17.247	6.993	31.473
beta0_H[13]	10.117	10.633	-10.784	10.186	30.226
beta0_H[14]	7.144	11.303	-15.479	6.974	30.492
beta0_H[15]	11.027	115.940	-215.147	10.502	253.224
beta0_H[16]	7.478	30.896	-57.211	7.980	70.946